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ABSTRACT 


The designer of stiffened cvlindrical pressure vessels 
is often faced with the problem of comparing structures 
of wiedly varying geometries. To do this, he makes use 
of a number of dimensionless parameters whose function is 


to reduce the number of variables in the design problem. 


The object of this thesis was to investigate the 
effect of the dimensionless parameter, А, on stiffened 
cylindrical pressure vessel design. In particular, the A 
region where colfapse of vessels due to yield-initiated 
buckling between the frames was identified and delimited. 
This region was used as the first iteration in a design 
procedure aimed at producing a vessel of maximum structural 
efficiency. Design curves which may be used to develop 
such a vessel are provided for shell thicknesses up to 
h/Dm 7» .01. These curves are based on both maximum 
principle stress and Hencky-von Mises stress criteria and 
provision is made in a computer program for use of other 
metnods. A modified Foppl Formula is used for frame 
evaluations. 


The formulas used range from conservative to liberal 
in their approximations of collapse pressures. However, 
use of even the most conservative formulas results in 
Structural efficiencies which compare favorably with 
present-day practices, 
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NOMENCLATURE 


s 4.72 
Cross sectional area of frame, іп 


Cross sectional area of frame plus plating of 
width b, in 


А1,А2,А; Plasticity coefficients, see Reference (10) 


b 


C 


Width of frame at shell, in 


Distance from neutral axis of bending to extreme 
outside fiber, in 


Outside diameter of shell, in 
Depth of rectangular frame, in 


Diameter to neutral axis of frame plus plating of 
width b, in 


Diameter to neutral axis of shell alone, in 
Modulus of elasticity, lbs/in' 
Eccentricity; in 

Secant modulus of elasticity, lbs/in' 

Width of flange of T frame, in 

Shell thickness, in 


Moment of inertia, frame plus plating of width b, 
S 
дү 


Moment of inertia, frame plus plating of effective 
width, see Reference (10), MEE 


Moment of inertia, frame alone, TN 

Moment of inertia, frame plus plating of width L, EM 
Unsupported span of plating, in 

Bulkhead spacing, in 

Frame spacing, in 


Unsupported span of plating taken as frame sapcing, in 





o 
o 





n 


Pc 


P1 


Number of complete waves in buckled configuration - 
clrcumferentially 


Collapse pressure, lbs/in? 


Actual collapse pressure desired, lbs/in* 


Pcli;Pc2;Pc3 Temporary values of po for iteration 


НЕ 


R 


Re 


Rm 


Ry 


V 


Wt 


OF 


Om 


Collapse load per unit axial length of a ring frame, 
lbs/in* 


Outside shell radius, in 
Radius to neutral axis of frame, in 
Radius to neutral axis of shell, in 


Radius to neutral axis of frame plus shell of width 
En 


Pome Cisplaced iy cylinder per unit length, iets 
Weight per unit length, lbs/ft 


Weight of salt water displaced by cylinder per unit 
length, lbs/ft 


Шек ency = defined in Section III 
Slenderness Ratio - defined in Section I 
Poissons ratio 

Axial stress at mid-thickness, lbs/in^* 


Tangential (circumferential) stress at mid-thickness, 
о un 


Material yield stress, lbs/in* 


Pressure factor - defined in Section I 





I. INTRODUCTION 


Stiffened cylindrical pressure vessels under external 
hydrostatic pressure fail by one or any combination of 
three collapse modes. While both theoretical and 
empirical expressions have been developed to determine 
collapse pressures for all modes, it is up to the designer 
to determine what structural propottions will result in 
the most efficient use of material strength and stiffness 
properties for a required collapse pressure. Two dimen- 
sionless parameters which have been widely used to define 
different geometries and their capabilities are w and 
A. V is called pressure factor and it is a ratio of a 
structure's actual collapse ure to the collapse 


pressure of an unstiffened tube with the same h/Dm ratio. 


Pc 
Y = 
h 
2---0 
D. Y 
A is called slenderness ratio. t is analogous to the 


slenderness ratio, l/r, of column theory, however, it 


includes terms which allow for different material properties. 





Despite ineir fairly wide use by designers in 


comparisons of different structural geometries, this 





researcher is unaware that any attempt has been made 

to develop a relationship between these parameters which 
would indicate their effect on structural efficiency. 

Also, it has been stated that these parameters are 
Sometimes misleading (7), however, no definite limits 

on their usefulness have been established in the litera- 
ture,and only a few other dimensionless parameters have 
been developed to be used instead. A need exists therefore, 
to define the limits of usefulness of both Y and À and 

to develop within those limits, relationships between 

them and structural efficiency which may be readily 

applied in the early design stages of a stiffened cylindrical 


shell. 


In this report, two failure criteria, maximum circum- 
ferential stress and Hencky-von Mises, have been used 
to develop relationships between À, y, and efficiency. 
In the course of the development, certain trends of the 
parameters have been noted and in the concluding remarks, 
useful limits of A and Y are indicated. In particular, 
the "low lambda" region has been defined and an intro- 
duction made to the concept of maximizing structural 
efficiency within that region. Design curves are presented 
to enable optimization within the low lambda region of 
structures with collapse pressures up to 2500 psi and examples 


of the process are given. Special note should be taken 





of Appendices B and C where remarks concerning very 


low lambda design are made. 





EDS EDOPMENT OF DESIGN CRITERIA 


Historically, the largest and most critically 
designed stiffened cylindrical pressure vessels have 
been intended for use as submersibles. Such vessels, 
under external hydrostatic pressure, have been known to 
collapse in one, or a combination of three modes. Since 
the mechanism which precipitates failure in each mode 
is different, the designer must take all three modes 
into consideration when designing for a given collapse 
pressure. However, he has a choice in the design method 
he may use when developing a highly efficient structure. 
He may design toward one mode of collapse, seeking to eli- 
minate the possibility of collapse in others while changing 
Structural proportions to obtain a high efficiency. His 
alternative is to simply design tovard the most efficient 
structure, not worrying about the mode of collapse as long 
as the required collapse pressure is obtained. 

An example of the first method may be found in the 
design of the relatively shallow-diving submarines of 
World War II. These vessels were designed to collapse by 
asymmetric buckling between frames. It was felt that 
imminent collapse vould be signalled by the appearance 
of lobes girdling the hull and the crew would be given a 


chance to surface before total collapse occurred (22). 





This explanation indicates that for the volume-limited 
submarines of that era, efficiency criteria could be 


established to a large degree by safety considerations. 


For deeper diving, weight-limited submarines, 
the first method proved inadequate and the second method 
was refined in a number of ways to allow development of 
highly efficient vessels. A refinement which will now 
be discussed is one which assumed that Y is limited to 
1.0. This assumption used the maximum circumferential 
stress criterion and minimized weight by using the 
smallest possible number of frames. It resulted in a 
structure proportioned such that À = .8. In order to 


visualize this approach, figure 1 has been developed. 


Figure l is a plot of Y versus d/h (frame depth/ 
shell thickness) for four A values. AS implied by 
the abscissa scale which indicates relative frame strength, 
the frames are rectangular and only their depth is varied 
to vary strength. The basic geometries used in this fig- 
ure are taken from a report by Trilling (19) and one of 
his curves (curve 2), as well as his experimental data 
(crosses) are plotted. The other three curves plotted 
on figure 1 represent the effect of changing A. In all 
three cases, 'he distance between frames was the only 


parameter varied. Values of Y were obtained by using 
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von Sanden and Gunther's formula #92A to predict collapse 


pressure. 


it will be noted that curve 3 of this figure lies 
almost entirely along the line, ф = 1.0. The A value 
of this curve is .814, and it is this relationship which 
forms the basis of the design method cited above. Just 
as the earlier method of designing for a particular 
collapse mode proved inadequate, this method has also 
been unable to develop the efficient structures ПН 
to resist high pressures. The remainder of this section 
will be devoted to developing a design method which will 


result in highly efficient structures. То begin this 


development, a further look at figure 1 is necessary. 


On the basis of curve 2 and the experimental points, 
one of Trilling's conclusions was that "except for very 
short thick tubes, . . . according to formula $92A the 
stronger the frame, the weaker the vessel." This was 
explained by noting that very strong frames form "hard 
spots" when, under external pressure, they contract less 
than the shell. The hard spots, in turn, cause stress 
concentrations in the otherwise uniformly contracting 


shell and lower collapse pressures. 


ENG 4 indicates, however, that this conclusion is 





not necessarily true, even for thin shells and large frames. 
The only difference between the model of curve 4 and 
Trilling's models is the decreased frame spacing which 
results in different A values; .898 for Tilling's models 
(curve 2) and .683 for curve 4. According to curve 4, 
shell strength actually increases with frame size when the 
frames are relatively close together. The seeming contra- 
diction of Trilling may be explained by an extension of 
his logic. If the frames causing stress concentrations 
are closer together (low lambda) the overall stress 
pattern will be uniform enough that premature yielding 


and buckling will not take place. 


The obvious lower limit of А is obtained when the 
frames are in contact (A-*0) (24) thereby creating a thicker 
shell (rectangular frames) or a webbed sandwich-type 
construction (T-frames). Both of these geometries result 
in putting the frame-caused stress concentrations in such 
close proximity that the overall stress pattern approaches 
uniformity. This uniformity of stress pattern implies that 
the yield strength of most of the material in the structure 
will be utilized before collapse occurs, which, in turn, 
implies that a structure of very high efficiency is possible 


лао À area. 


With a lower limit of À established at 0, an upper 





limit is necessary to define the range of geometries 
where the effect of increasing frame strength is not 
detrimental to collapse pressure. This limit will be 
meterred to as Am. For structures with À values above 
BEncereasing the frame strength weakens the vessel. 
For A values below RT increasing the frame strength 
strenthens the vessel. An initial value of this limit 
can be found by direct use of Windenburg's approximation 
of von Mises' formula for collapse between frames due 
to elastic shell instability. The development of this 
formula is outlined in Appendix A. In its simplest 


form, this formula is: 





ф = De 
‚2 
Recalling that у = Pc , Windenburg's formula 
h 
E 


determines instability collapse in terms of the hoop stress 


relation. i.e., 


hc Las 
Po N; “p. у 2 
therfore, when 1.3 = 1.0, failure is due primarily 


x 2 
to yielding of the material rather than plastic instability. 


u iS point, A = 1.14. If A 1.14, then instability 


collapse is indicated since Y 1.0. In this situation, 





-10- 


strengthéning the frames cannot delay collapse. It may 

be tentatively assumed that when A l.14, strengthening 

the frames increases overall strength. If this is so, 

then an initial approximation of Ав as 1.14, TË is important 
to note that this is a very approximate method for 

Are mining An. A glance of Trilling's results on figure 1 
will indicate that collapse at Y values below 1.0 has occurred 
in structures with A values as low as .89. Also, with 

this method, Am 1s independent of any individual structure 
parameters. Since this is questionable, another means 


of determining Am is necessary. 


Returning to curve 3 on figure 1l, it will be recalled 
that for the A value used in that curve, ф = 1.0 for all 
frame sizes. As long as the frames are of adequate 
Strength or stiffness to carry a given percentage of the 
shell load, their size has no bearing on the strength of 
the cylinder. Figure 2 elaborates on this point. The 
curves of figure 2 are collapse pressures, determined 
by 592A, plotted as Y against values of A. Each curve 
represents one size of frame and all structures have the 
Same h/D' ratio, so that only the distance between frames 
has been varied to obtain changing À values. The point 
at which all curves cross is at ý= 1.0 and represents 


the afore-mentioned Ал. Unlike the approximate Windenburg 


determination of this point, in this determination, Ag 
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depends on the value of h/D as well as the fact that 
ф = 1. This dependence is borne out in the derivation 
of an expression for Am from formula #92A. This deri- 


vation is found in Appendix E and results in the 


equation: 
D = 2.055D eN (A) 
m man a) A 
h E 


which was used to plot h/D, versus A, for different 


values of бу on figures; 


Figure 3 may be used to determine for a given 
Dp and O what A is necessary to prevent collapse be- 
tween the frames due to instability (y < 1.0). Areas 
to the right of the curve in question are areas of 
instability collapse, those to the left, primarily 
yield failure. As an example of this, Trilling's 
experimental results have again been plotted. Due to 


slightly varying geometries and o all the models he 


y $ 
used fall within the areas of the dots shown. The 
figure indicates that for с = 30,000 psi (see appendix) 
all the models failed at a i$ below 1.0. Reference to 
figure 1 shows that this was indeed the case,and a 
further glance at Reference (19) indicates that a number 


of Trilling's models failed by lobar buckling or general 


instability, rather than by yield initiated collapse. 
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In Appendix A, a number of formulas and yield 
criterion are cited as being more accurate than the 
maximum circumferential stress criterion in determining 
collapse pressure due to yielding. One of these, Wenk's 
formula, has been investigated using the same models 
used in figures 1-3 in order to ascertain the influence 
of increased accuracy on the computationof i,. Figure 


4 is the result of this investigation. 


Wenk applied formulas £92 and #92A to the Hencky- 
von Mises yield criterion at mid-thickness, mid-bay, 
and obtained collapse pressures within nine percent of 
actual experimental values. The Hencky-von Mises criterion 


in two dimensions is: 


2 2 2 2 
су = QOmq^ - OgpOp + Op 


Bm at is applied to an unstiffened shell where Om = 

20p, yield failure occurs at Om = 1.160y+ i.e., Ü = 

1.16. Just as figure 2 indicated that for the maximum 
circumferential stress criterion, there is a À value 

where frame size has no effect on cylinder strength, figure 
4 shows that the same is true for the Hencky-von Mises 
criterion. However, the À value in this case occurs when 


I y6 instead of 1.0. MWenk's formula: 
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\r: cao cS 
like formula #92A, is actually the hoop stress relation 
times a factor(defined in Appendix E), When it is solved 
for F by setting Y = 1.16, a An versus h/D,, relation 
Similar to formula (A) may be found. This formula is 
shown below and a detailed derivation of it is given in 


Appendix E: 
(B) 





If plotted in the form of figure 3, this equation would 
result in a shift of all curves .11A to the right of 

those obtained from #92A. This shift reflects the conser- 
vative nature of von Sanden and Gunther's formula 


and criterion. 


Thus far, equations and criteria have been developed 
to determine the extreme values of A between which 
stiffened cylindrical shell failure is primarily due to 
yielding of the material between the frames. Justification 
for using this as the design collapse mode comes from a 
number of sources.  Wenk, a major advocate of the "one- 
hoss-shay" collapse philosophy, states that the approach 
of designing for yield failure is not inconsistent with 


that philosophy (221. Because of the effects of such 





E 


unknowns as material residual stresses and out of roundness, 
design for collapse by all three modes simultaneously 
ШЕГЕШЕТІСІЛЕ, if not practically impossible. Therefore, 
the mode which is least affected by such strength modifiers 
should be the design collapse mode. Conversely, the 
collapse mode most affected by the modifiers should be the 
first mode taken into consideration and designed out of 
existance. Out of roundness is the primary reason for 
local failure of frames which can precipitate premature 
general instability collapse (15). For this reason, 
collapse by general instability has been noted as being 


of prime concern in the preliminary design (10). 


Reference to formulas for general instability given 
in the appendix will reveal that this collapse mode, unlike 
the others, is a function of cylinder length (between 
bulkheads), as well as frame area, moment of inertia, and 
Spacing. Judicious selection of bulkhead spacing (if a 
free variable) and minimum frame sizes can eliminate 
this collapse mode from consideration. This should be 
done in the first stages of design development. In the 
final stages collapse pressure for this mode should again 


be checked and out of roundness taken into consideration. 


With general instability designed for, collapse 


between the frames due to buckling should be eliminated 





-18- 


Бу insuring that the final design À will be between 0 

and Am. This means that the yield strength of the material 
will be utilized as fully as possible by insuring against 
premature failure initiated by instability. In this way, 
the benefits of using high strength materials will not 


be nullified by instability collapse at low pressures. 


In the next section, considerations for the final 
step of designing for yield initiated failure between the 


frames will be developed. First, however, the question 


of frame adequacy must be looked into. 


Figures 1, 2, 3, and 4 were all developed from 
analysis of cylindrical shells with rectangular stiffeners. 
T-shaped stiffeners are in general use today however, so 
the question may be raised as to what effect the shape 
of the stiffeners has. An investigation of the formulas 
for yield failure previously cited will reveal that they 
are all based on the assumption of adequate frame strength. 
The only frame parameters émployed in them are flange 
width in contact with the shell and’ frame area. At А = 
КОШ и песе parameters cease to influence collapse 
pressure. Therefore, it can be stated that frame shape 
has no influence on formulas (A) or (B). In general, 
as long as adcquate strength is maintained, frame shape 


does not influence the collapse pressures of stiffened 





Ec 
cylindrical shells. 


The question of what constitutes adequate frame 
strength has been directly approached by a number of 
investigators (2,11,19). The standard practice of 
von Sanden and Gunther (20) and other early designers 
was to design the frame to be able to withstand 1.1 
times the total pressure on the length of plating it 
supported. This procedure, which employed the Foppl 
formula, was found to result in "excessively heavy frames" 
(19). In addition, the Foppl formula has been determined 
to give conservative results for frame buckling pressures. 
In light of these facts, the frame design requirement 
chosen for the purposes of this Desert was that the frame 
be strong enough to support no more than the load on 
one effective length of shell. The conservativeness of 
the Foppl formula was accounted for by assuming that it 


predicted frames alternatively 4, 5, or 6 times too 


strong. The resulting formula used is: 


24EI, 
DE rame е рос әлегі = — X (4,5,0r 6) 
Dy Ly 


Before proceeding into the actual design iteration 
resulting from the establishment of design criteria for 
stiffened cyl.ndrical pressure vessels under hydrostatic 


pressure, a recapitulation of the criteria is in order: 





=2:0= 


Upon definition of -operating requirements, initial 
Structure parameters should be established to 
minimize the possibility of general instability 


collapse of frame and shell. 


Possibility of collapse due to instability 
between frames should be reduced by limiting 
design choices to geometries located between 


А = 0 and A = Ат 





oue 
ШЕРУ VESSED EFFICIENCY 


The measure of efficiency used in this paper 
was suggested by Wenk in (22). It is slightly more 
complicated than the traditional Weight/Displacement 
ratio in that it takes into account collapse pressure. 
PoV 
10%w 
р. is the collapse pressure - lbs/in2 


V is the volume per unit length of cylinder - in*/ft length 


W 1s weight per unit length of cylinder Іре/ге length 


When used in a qualitative manner, this term clearly 
shows the relative efficiencies of different geometries. 
However, it is not independent of material properties 
such as density and yield stress. The efficiencies in 
this report have been calculated using HY 80 as the 
material. In order to apply them to materials of other 
densities and yield strengths, they must be multiplied 


by the folloving factors: 


O . 284 
30, 000 and/or — ES 


For example, a geometry vhich would have an efficiency 
of 1.25 for HY 80, has an efficiency of: 


50,000 ,.284 _ 94 ig 


D aaoo 
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when the material is 50,000 psi aluminum. Once found іп 


the above manner, the efficiency may be used to calculate 
the traditional Wt/Displacement ratio using the following 


conversion formula: 


Wt Ре 


À 3720n 


Reasons for this manipulation of n will become more clear 


at the end of the next section. 





De 
IV. DESIGN ITERATION 


In the preceding section, design criteria have 
been developed which, if used, will produce a generalized 
cylinder geometry whose Avalue will be between 0 and Хр: 
This will result in a collapse pressure factor in excess 
of 1.0, thus insuring that collapse will be primarily 
due to yield and the yield ME of the material will 
be highly utilized. This section will outline development 
of an iterative method to find the particular А апа vy 
values which result in maximum structural efficiency 


for a given collapse pressure. 


Ficures 5, 6, and 7 are plots of efficiency against 
E terent values of h/D. and frame strength factors. 
As stated previously, Foppl's frame formula has been used 
in a modified form vhich arbitrarily reduces its apparent 
conservativeness. In determining frame size for figure 
5, Foppl's formula was multiplied by a factor of 4. In 
figure 6, the frame factor was 5 and in figure 7, it was 
6. The effect of this increasing theoretical strength 
of frames was to reduce frame size and increase apparent 
structure efficiency in succeeding figures. Since this 
is an arbitrary change, model tests are necessary to 
determine the actual frame factor which should be used. 


It was felt that a fange of factors such as that covered 
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would indicate design trends which may be beneficial in 
pEeucersl-cuseswsthout the necessity of establishing rigorous 
theoretical frame design criteria. Further substantiation 


of this method will be provided later in this section. 


The formula used to determine collapse pressures 
for these three figures was r92A (maximum circumferential' 
Stress criterion). Р tHis is now considered to 
be less accurate than a number of possible formulas, it 
has the advantages of always being conservative and 
having been widely used. The curves in figures 5, 6, and 
би всаће that for a given h/D, and frame factor, the 
efficiency increases with Auntil a maximum, Nmax’ 15S 
reached. Thereafter, the efficiency steadily decreases. 
with increasing А . It is this point, Mmax» that the 
designer of stiffened cylindrical pressure vessels is 


interested in, therefore, only enough points have been 


plotted to indicate the hump of each efficiency curve. 


An interesting trend which is readily apparent in 
all three figures is that the efficiencies of all h/D 
combinations lie on the same line at the point where 
A = Am for each h/D (see figure 3). In figure 5, this 
line has a slope of about -.5. In subsequent figures, the 
slopes appear the same, so, since the lines begin at higher 


values, the final efficiencies will be higher. Of course, 





ae 


extension of the lines upward to maximum efficiency also 
requires extension to the left to A - 0 and extension of 
the series of h/D curves to extremely thick cylinders. 
This extrapolation of computed data can be justified on 
the grounds that only very short, thick cylinders fail 

by pure yielding. Failure in this way takes the fullest 
possible advantage of material yield strength and results 
in the highest possible structural efficiency. Since 
efficiency in this report is also a function of displace- 
ment to weight ratio, the veracity of that PS 


Statement, in the general sense, is dependent on material 


density. 


Another point which is illustrated in figures 5-7, 
ашпас гог еасһ h/D,, there is only one value of À at 
which maximum efficincy is reached, regardless of frame 
factor. The effect of the frame factor 1s to change the 
actual maximum efficiency reached, but it has no effect 
on the A value where it is reached. This direct relation- 
ship between h/D n and A at maximum efficiency is shown 
190ге 8. It may be explained in part by the fact that 
K I A nor h/D, contain any frame size parameters and 
both Wenk's formula and formula 792A assume adequate 
frame strength. Also at that point of maximum efficiency, 
there are dis.rete values of W which are dependent on 


frame factor. In order to show the relationship of frame 
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factor, A and Y at maximum efficiency, figures 9 апа 


10 have been plotted. 


Before entering into the use of these figures, 
however, a slight digression is necessary.  Recalling 
that formula 192A, used to obtain figures 5-7, was 
shown to yield less accurate collapse pressures than 
other formulas, it became necessary to investigate a 
formula of higher accuracy to determine if the same 
trends still held. The formula used was developed by Wenk 
ena has been introduced in Section II. Data similar 
to that shown in figures 2, 6, and 7, was developed for 
the same geometries using Wenk's formula. The actual 
results are given in Appendix G. Integrated curves from 


these results have been drawn on figures 8 and 10. 


Figure 8 is a graph of h/D versus A for points of 
maximum efficiency. The two curves represent results 
of the use of formula 292A and Wenk's formula. It will 
be noted that Wenk's formula shifts the curve to the right. 
Since this expression has been known to predict collapse 
pressures in excess of actual values, it may be expected 
that the curve for actual structures will lie between the 
curves shown. No mention of frame size is made on this 
graph, as in neither formula did frame strength alter 


the relationship of h/D and A at: maximum efficiency. 
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Figures 9 and 10 are graphs of Mg ж апа Y Lor 
:92A and Wenk's formula respectively. The three lines 
in each case indicate the three frame factors used. 
Comparison of these figures shows that increased accuracy 
does not change the trends previously noted. In fact, 
the n versus À curves are exactly alike. The 25 percent 
increase in efficiencies in figure 10 over figure 9 is 
another reflection of the conservativeness of formula 


192A. 


With development of figures 8-10 complete, an 
iterative method to determine geometries of maximum 
efficiency using these figures can be outlined. Geometric 
parameters which can be found are h/Dm and À. Given 
material properties, distance between stiffeners may 
be computed. Actual stiffener parameters are left up to 
the discretion of the designer as the only requirement of 
the formulas used in this procedure is that the stiffeners 
be of adequate strength and stiffness. For comparison 
purposes, the frame parameters used in this report and 
those frame dimensional ratios coinciding with maximum 
efficiency structures are given in Appendix C. To begin 
the iterative procedure, the information vhich must be 


specified is: 


1. type of analysis to be used 


a. maximum principle stress, #92A 





Em 


b.  Hencky-von Mises stress at mid-thickness, 
mid-bay, #92 and 492A 
Ри" тате factor to be used: 4, -5, Or 6 
3. ultimate collapse pressure desired 
4. material properties, бу апа Е 
(боли ресе о шила E == EX LON) 


5, outside diameter 


— m 


An example of the procedure is given on the next page. 
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Example 1: 


Given: Use of 492A, Frame Factor - 5 
Sy = 80,000 psi 
D ee Obert. 
p, = 1000 psi 

oteps: 


ПА слап цене ош тоорио ево 





pe = 2(h/D)o, h/D = h/Dnm 
(h/Dm) = 1000/ (2*80,000) 
(YD 29006 92 
À = .48 from figure 10 
V = 1,29 from figure 11 
y = Pc 
= 1.29 
е 


Dene 00290" psi which as too high 


2. Reduce h/Dm by division 


Pci = 1000 = 775 psi 
E229 
(h/Dm) = 775/(2x80,000) 
(h/D, ) = .00484 
À = .58 
Y = 1.41 
Po = 1.41*775 = 1090 psi - which is 


still too high 





Answer: 


Results: 


36-2 


Reduce h/D, again 





Рс? 7 oom an PS 

(h/D,,) = 710/(2x80,000) = .00443 

A = ,.606 

ф = 1.44 

Ре = 1.44x710 = 1020 psi - which 


is close enough, but try a 4th 


iteration 
Reduce h/D, 
710 
Рез 7 31.02 7 9?" Psi 
(h/D,,) = 697/(2x80,000) = .00436 
À = „613 
Y MaS 
Ре = 1.45*697 = 1012 psi 
À | = с 
h/D = 00436 
n = 1.22 
Y а 5 
р = ЛО FE. 
h 2x xin: 
Le "БЕГЕ. 
Wt. MTS 000 БУТЕ 


conventional Wt/Disp = .222 
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foe e toCUSSION OF RESULTS 


The results of this report may be divided into two 
classes; those due to analysis and those due to synthesis. 
Analysis by formula #92A and Wenk's Formula at low 
lambda values has shown that considerable gains in 
efficiency can be realized by placing frames close together 
and reducing shell thickness. This results in | values 
greater than 1 and helps to realize the stress potentials 
of high strength materials. The extrapolation of the 
design curves in figures 11 and 12 appears to justify 
the use of sandwich construction for very high pressures. 
Such construction is analogous to design in the very low 
lambda region, and a discussion of this is included in the 
appendix.  Inferences which may be drawn from the treat- 
ment of very lov lambda design in Appendix B concern the 
apparent inapplicability of A and bin this region. This 
results from the fact that in this region, frame strength 
contributes so much to overall strength that it cannot 


be ignored as it is in both A and ү . 


At the high end of the low lambda region, maximum 
structural efficiency can be realized only in relatively 
low pressure vessels. This is shown in two examples in 


the appendix. 





Eogr 


The result of synthesis in this report is the 
iterative method outlined in Section IV. By its very 
nature, this method results in thin-skinned vessels 
with frames relatively close together. The questions 
which this geometry may raise concerning modes of collapse 
other than yield-initiated, must be answered by indepen- 
Gent analyses. It will be noted, however, that the 
geometries used to plot the curves of figures 8-13 were 
analyzed for failure both by instability between frames 
and general instability and the indicated collapse mode 


for all was yield-initiated failure between rings. 





0 = 


ED CONCLUSIONS 


1. The dimensionless parameter, A, may be used to 
define all stiffened cylindrical pressure vessel geometries 
except those in which the frames provide little support 
because of their wide spacing, and those vessels in which 
the frames provide most of the strength by virtue of their 
size and close spacing. " 

2. Absolute structure efficiencies in the low lambda 
region are highly dependent on frame spacing and only 
slightly dependent on frame size (as long as the frames 
are of adequate strength). 

3. Relative structure efficiencies in the low lambda 
region are virtually independent of frame size. The lambda 
values at which maximum efficiency is reached depend on 
A DO ratio and are not affected by changes in frame size. 

4. Structures with w values well in excess of 1.0 
are reasonable and generally provide increased efficiencies 
over those where y is limited to 1.0. 

5. As the design collapse pressure increases, the 
value of lambda which coincides with maximum structural 
efficiency decreases. This implies that higher efficiencies 
can be obtained by reducing shell thickness well below that 
required by the hoop stress relation and simultaneously 
increasing the number of frames. 

6. At extremely low lambda values, general instability 


becomes the most likely mode of failure. 





E а 
VII. RECOMMENDATIONS 


Since this thesis was limited to studying the 
parameter, А, and its effect on efficiency, other 
definitive parameters, such as those now used at the 
David Taylor Model Basin should be investigated for 
their effect on efficiency. It has been noted that the 
lambda term is relatively useless when the frames are 
very close together and are very strong. This condition 
Should be corrected by defining a new parameter similar 
to lambda in its application, but allowing for the effect 


of frame strength in geometric comparisons. 


The basic approach of this thesis should be used 
to extend the curves of maximum efficiencies, found in 
Section IV, into the regions of higher collpase pressures. 
With this extention, a thorough re-examination of general 
instability collapse in these regions should be undertaken. 
Finally, experimental tests of maximum efficiency geo- 
metries, as developed by the design iteration, should be 


undertaken to provide evidence of the method's worth. 
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APPENDIX A 


HISTORICAL BACKGROUND 


Prior to 1917, the design of stiffened cylindrical 
pressure vessels was based on the assumption that vessel 
strength depended primarily on frame strength. The shell 
was treated as a flange of the frame and strength was 
calculated from equations which defined collapse by 
instability of frame rings under radial toads. A widely 
used formula for this was developed by Levy, but is 


commonly called Foppl's Formula (19,24): 


While this formula was relatively accurate for the 
time, and where frames alone were involved, when it was 
applied to a stiffened cylindrical shell, its shortcomings 


became apparent as higher collapse pressures were required. 


In order to improve upon this frame-oriented theory, 
von Sanden and Gunther derived expressions based on the 
m@ieservation that stiffened cylindrical shells were not 
Ferlleetıons of connected rings, but, in fact, were closed 
e linders reinforced by ring frames. Deciding that the 
shell was the real strength member, they predicted failure 


due to shell yielding at the area of maximum stress. This 





SASE 


maximum principle stress criterion could be satisfied 
by the longitudinal stress at the frames or by the 
circumferential stress midvay between them. Formulas 
for both stresses were developed and the well-known 
formulas #92 and #92A resulted when pressures were 

for in terms of the stresses. Substituting material 
yleld stresses into #92 and #92A resulted in two yield 
pressures for the structure, m the lower one was 
said to predict failure. It must be emphasized that 


failure was predicted on the basis of the onset of 


yielding. 


For a number of years after the 1918 publication of 
these formulas, $92 was thought to give the critical 
collapse pressure since the greater stress was longitudinal. 
In 1935, however, this assumption was found to be in 
error when circumferential stress was determined to be 
critical in the collapse of thin-shelled pressure vessels 


(19). 


The von Sanden and Gunther formulas are: 


_ ZUG 
Py = e е ВЕНЫБЕБИСИНЕЭРНИСТ 
+ 1.81К [.85-вВ 
ЕБ 


#92 





492A PE 2h /p Су 
IN ср 
1 + 6 


The variables K, H, B, and B are defined in the Formula 


section, 


About this same time, it was generally recognized 
that the area between stiffeners very rarely collapsed 
as a result of pure yielding. In anything but thick 
cylinders, collapse occurred from a combination of 
yielding and buckling when the distance between stiffeners 
(effective length) was below a critical length. This 
critical length was dependent on shell L/D and h/D ratios 
as well as material properties. When the distance between 
Stiffeners was above the critical length, the tube was 
considered to be of infinite length and the formula of 
Breese and Bryan for collapse by instability could be 


applied: 


5 = 2E 
c ` m Lan 3 
——À3 (by) 


For cylinders below the critical length, formulas 
for collapse by instability were developed by von Mises 


and Southwell. Windenburg simplified von Mises lengthy 





Aly 


formula and, by sacrificing.a very small degree of accuracy, 
came up with the Eollowing formula (25): 


DID (ШОШ Б 
p — tn mn 


м 500601 
At the time Windenburg developed the above formula, 
a Pe versus L/D coordinate system was commonly used to 
plot experimental and theoretical results for changing 
values of h/D. In an effort to develop an analog to the 
Slenderness ratio, L/r, in column theory, and ME SMOD 


of results, he developed a new coordinate system based 


on the coordinates: 


ur р 

р, = и 
2h/Dm 

de = | (0)? 


(gon Da) ° 


Windenburg then corrected his coordinates for the physical 
properties of materials in the following manner and came 
up with the well-knovn variabels, slenderness ratio А, 


and pressure factor, V . 
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This done, Windenburg rewrote his instability formula in 


terms of his new variables: 


O 
№ = 
A*-e 
ЗА 
_ 2 
where e - .045 19996, E =) 4 
Е 91 100t/D 


may be neglected in general to get 


1539 
Y = 


у 2 


This formula is plotted on figure À-1 and may be 
compared there with the plot of von Mises' formula and 
numerous experimental results. One thing which may be 
seen at this time is the good agreement these formulas 
have with eachother and the experimental results above 


a A value of 1.4. 


After a decode of relative inactivity, the investi- 
gation of stiffened cylindrical pressure vessels was renewed 
in the early 1950's. The main intent of the investigations 
of this period was to correct the number of theoretical 
deficiencies in formulas +92 and #92A. The first two 


deficiencies were that the formulas predicted linear 
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relationships between pressures and stresses and ignored 
the effect of moments developed due to pressure at the 
ends of closed cylinders. Both problems vere corrected 
by Salerno and Pulos (17), who developed stress formulas 
similar to #92 and %92А but taking both the above 

moments and the influence of axial stress in the shell 
tending to expand the frame into account. Their analysis 
resulted in a non-linear ship between pressures 
and stresses which, however, did not alter the results 


of stress calculations significantly from those found by 


#92 and #92A. 


The next problem found with von Sanden and Gunther's 
theories had to do with the inadequacy of the maximum 
principle stress criterion. Use of the Hencky-von Mises. 
preld criterion: 


2 ~ 2 2 
O = O = Ото +0 
y T IE F 


Be exterior mid-bay point resulted an predictions closer 
to actual collapse pressure than obtained with the former 
criterion. (9). However, here again collapse pressures 

were associated with the onset of yielding and no allowance 


for plastic reserve strength was made. 


This condition was corrected by Wenk et. al. (21) in 


a somewhat haphazard manner by taking the Hencky-von Mises 





= 


yield criterion at mid-bay and mid-thickness. Although 
Wenk termed this a "necessary, if not sufficient empirical 
yield condition to produce failure", experimental results 
have since indicated that thought this formula predicts 
within nine percent of actual collapse pressures, it 

tends to predict high. (9). An interesting observation 

on the use of the Hencky-von Mises criterion is that for 
an unstiffened cylindrical shell, от = 20p + This results 
in failure of such a shell occurring at oq - l.160y; 


stated in Windenburg's variables,v = 1.16. 


The use of formulas $92 and 292A in the development 
of Wenk's formula appears to be the last time these 
long-standing elastic stress determinators were used in 
a yield criterion. The elastic stress formulas of 
Salerno and Pulos (17), having been recognized as more 
rigorous and slightly more accurate, were used in the next 
collapse criterion developed by Lunchick.  Lunchick 
utilized the Hencky-von Mises yield criterion, however, 
he specified that collapse did not occur until the load- 
carrying capacity of the shell vas completely exhausted 
at the critical mid-bay point.  Lunchick's three-hinge 
mechanism stemmed from his assumption of plastic hinges 
forming at the frame at the same time a completely plastic 
hinge had developed at mid-bay.  Lunchick's method predicted 


collapse pressures with better accuracy than any other 
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in the range of low values. However, its inclusion of 
design curves to determine cylinder plastic reserve 
strength has made it more difficult to apply in computer- 


based design methods. 


It is now generally accepted that stiffened cylindrical 
Shells fail by one of three distinct modes or any 


f 


combination of them. These modes are (14): 


Mode 1. Axisymmetric collapse of the shell between 
adjacent ring frames. 

Mode 2. A a collapse of the shell between 
adjacent ring frames (lobar). 

Mode 3. Overall asymmetric collapse of the shell 


and frames together. 


Mode one had traditionally been referred to as 
axisymmetric shell yielding.(9,22,23). It has been 
pointed out, however, that this mode is the result of a 
combination of yielding and axisymmetric buckling (14,22) 
and the main difference between this mode and mode two is 
the difference in resulting collapse patterns, not basic 
collapse mechanisms.  Axisymmetric collapse pressures 
(е been calculated with good accuracy through use of 
yield criterion alone though, and it is recognized that 


pre-collapse stresses, as determined by von Sanden and 
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Günter (20) or Salerno and Pulos (17), are the primary 


canse ОТ failure in this mode. 


The problem of axisymmetric collapse by elastic 
buckling has been solved by Salerno and Pulos in (13). 
Lunchick (10) has investigated axisymmetric failure by 


inelastic buckling and derived the following expression 





fOr it: 
2 A125? 
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von Mises' and Windenburg's formulas for the calcula- 
tion of failure pressure by elastic asymmetric (lobar) 
buckling between ring frames have been given. A more recent 
solution to this problem has been found by Reynolds and 
is outlined in reference (14). Reynolds has also developed 
a formula for the inelastic buckling of stiffened cylindri- 
cal shells which appears to be reasonably accurate, but 
eS Consistently high in its calculation of collapse pressures 


(1). 


Mode thr e, collapse by general instability, has 


been found to be of primary improtance in the preliminary 
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design of pressure vessels (14). Solution of elastic 
general instability collapse was first undertaken by 
Tokugawa, who was also the first to identify this mode. 
His formula for collapse pressure, (which was later 


confirmed by Bryant), is as follows: 


~ 


Eh ne 2 








_ EIe 
Pcr = En > (и D 
n2 — 1 + > ) n? + а?) 2 R’Lr 
TR 
where Л p 


The original Tokugawa formula as given in (19) has a 
number of extra terms in ıt which are seen to be 


negligible for geometries of general interest. 


The lengthy formulas developed by Kendrick for 
elastic general instability are given in (14), however, 
a graphical solution of Kendrick's equation as developed 
by Reynolds, and refined by Ball, is included as figures 
A-2 and A-3. Іп e solution, the sum of two effects, 
each represented by a graph, is used to determine collapse 
pressure. Figure 2 is the shell effect and figure 3 
the frame effect which reflects the bending stiffness of 
the cylinder. Figures 2 and 3 also indicate the Bryant- 
Tokugawa solution. Collapse pressure by either set of 


curves is found in the following manner: 
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In order to round out this development of stiffened 
cylindrical pressure vessel collapse theories, mention 
must be made of the effect of actual physical manufacturing 
conditions. In all but the smallest and most painstakingly 
fabricated cylinders, imperfect production processes are 
the cause of a number of collapse strength modifiers. 

Those modifiers introduced into the material itself are 

in the form of stresses due to fabrication, rolling, and 
welding of flat plating. In thick homogeneous plating, 
the homogeneity itself may be questionable and thus cause 
stress concentrations. A major strength modifier is the 
Euperfect circularity of most welded cylinders. This 
modifier has received some theoretical consideration as 

it has a considerable Are on the resistance of a stiff- 
ened cylinder to failure by general instability. A 


formula for this is due in part to Kendrick (2): 





IE 
PIR l+ B Ec Pl 
И. тотса с 2-е \ 
2 = 
ZB 1 + B Д PcrP] 


EN ong as theo... is below material yield stress, 


collapse by general instability is not indicated. 
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This section has been included to give a brief 
background of the development of theoretical collapse 
predictors for stiffened cylindrical shells. In particu- 
lar, a history of the variables of à and Ņ is helpful 


in understanding their use in this paper. 
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APPENDIX B 


DISCUSSION OF VERY LOW LAMBDA DESIGN 


As noted in Section IV, overall structural efficiency 
approaches its highest value as Y goes to 1.0 and A goes 
to 0. For T-shaped frames of the geometry used in this 
report, À values in the area below .2 may be considered 
very low lambda design. For other frame geometries, this 
may be higher. For example, the sandwich construction 
in figure A-4 has a À value of about .337. The frames used 
there have thicker webs than flanges, whereas the reverse 
ms true for the M mes of this report. The membrane shell 
design, and other complications resulting from the use of 
high strength titanium also plays a part in deciding the 


optimum À value. 


It was noted in Section IV that the danger of general 
instability collapse would be greatly increased by the 
use of a very low lambda design. This has been shown to 
be true in model tests of cylindrical shells with 
geometries similar to that of figure A-4. In fact, in the 
testing of these models, a new mode of collapse was observed 
- that of plastic general instability - and new theories 
had to be developed to predict it. It may also be observed 
that collapse of these models took place at Y value in 


excess of 3.0. This suggests that at very low lambda 





EG 


values, 6 , becomes a rather meaningless term unless 

the framing is very small in comparison to the shell 
thickness. Indeed, the large effect of framing at these 
A values throws into doubt the meaningfulness of À itself. 
This doubt has been substantiated but not fully explained 


in reference (6). 
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APPENDIX C 


FRAME PARAMETERS 


A T-shaped frame was used in the calculations 
leading to figures 5-10. The proportions of this 
frame were taken from a work by McGinley (11) and are 


as follows: 


~. 


frame width/frame depth = ]1.0 
frame width/web thickness = 8,5 
flange thickness/web thickness = 1.7 


Figure A-5 is a graph of b/D versus X for the frames 
for both formulas 292A and Wenk's., Since the geometry 
of optimized vessels was found to be a function of required 
collapse pressure,and not required diameter, the optimized 
frame dimensions should be proportional to the major 
cylinder dimensions, hence the b/D parameter. Figure A-4 
may be used in conjunction with figures 8, 9, and 10 to 
determine the frame dimensions. Care should be taken, 
however, to match yield formulas and frame factors when 


using these figures together. 
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EXAMPLES 


APPENDIX D 


OF DESIGN ITERATION 
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Example 2: 
This example is included to indicate the results 
of the iterative method when applied to relatively low 


pressure vessels. 


Given: Use of #92A Frame Factor 5 
бу = 30,000 psi 
Ре = 400 psi 
D = oO Eb. 

Steps: 


1. Initial h/D, guess from hoop stress 


(N/D) = 400/ (2х30,000) 

(h7D,,) = .00666 

À = .495 

V = 1.30 

Pc = 1.30*400 = 520 psi 


which is too high 


2, Reduce h/D, by division 





5:24 = 400 = 308 psi 
J: 
(h/D,) = .00512 
À = .565 
Y = 1:59 
Pc = 1.39*308 = 428 psi 
which is still too 
high 
Vact ш _428_ = 1.07 











Answer: 


Results: 


Enc 


Reduce h/D, agaln 





288 psi 


406 psi 


which is close enough 





Pc2 = 308 E 
= 1.07 
(h/D,) = .00480 
À = .584 
y = 1.41 
pr = 1.41%288 = 
A = 207 
h/D. = ,00480 
n ME соо 
80,000 
n = ,476 
Y = 1,41 
h =e, me 
n БИЕ 
Wt = Pe 
A 
37207 
Wt - 


PS 123 


406 
3720(.476) 





= 57 


Example 3: 
This example is exactly like example 2 except that 


a higher strength material is used. 


Given: Use of #92A Frame Factor 5 
бу = 50,000 psi 
Pc = 400 psi 
D = 40 ft. 
Answer: À = .825 
pop MEL 200052 
n = 1l.11* 50,000 _ 
80,000 
n = .694 
Y = 1.735 
Results: D = 40 ft. 
h АИ КО rn. 
L: = 1.83 ft. 
Wt = .157 
А 
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Example 4: 
This example is included to indicate the use of the 


frame parameter curves in Appendix C. 


Given: Use of #92A Frame Factor 6 
Oy, = 100,000 psi 
Ре = 1500 co 

Steps: 


то за  пур guess 


(h/D,) = 1500/(2*100,000) 

(h/Dm) ~ = .007 

\ = .48 

Ü = 1.33 

Dc = 1.33%1500 = 2000 psi 


2. Reduce h/Dgq 





Pcl = = — 1120 

(пир) = 1130/(2%100,000) 

(h/Dm) = 1.0056 

Х -z god 

Y O 

Pc SA A psi 
1580 ИЕ 


Vact 7 1500 





mIs 


ee Da again 


Pc2 
(H/Drp) 
(h/Dq) 

À 

Y 


Pc 


Answer: X 


(h/D,) 


= » 


Results: 


о O Jie 


E 


Frame Dimensions: 


b 


80,000 


Te - 1075 psi 
20757122100,000) 

. 00557 

тоо 

1.41 
1.41%1075 = 1515 psi 


close enough' 


MOS 


200537 
100,000 


4471 

.00248 
40 ft. 
2460 Ln. 


J.42 ft, 


2502 


ie om yn. 
in. 
ола. 


2705 5n. 


je 12 
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APPENDIX E 


FORMULAS AND DERIVATIONS 


In this appendix, the formulas used in the er 
programs in Appendix F are listed. Also, the derivation 
of formulas (A) and (B), Section II, and the solution 
of the Hencky-von Mises criterion for unstiffened tubes 


=, 


are presented, 


FRAME FORMULAS: The formulas below are used to 
predict frame failure by instability. Collapse of the 
frame by yielding was not investigated as the frames 
were not required to hold up after shell failure and the 
maximum stress in the frame is always less than the 
maximum stress in the shell (19). However, this stress 
relationship may not hold in the very low lambda region. 


Foppl Formula: 
24ЕТЬ 
Ре = ° 
Ор Ly 


Formula #88 
von Sanden and Gunther: 


24EI, (148) 
Pa = eee 





Du, °b bet .85В ) 


B 
bh 


where B z —————— 
А + bh 





ES 





L/Dm 
0 = 1.82 h/D,. 
| 4L NB 
B = 25 8 
N = coshd8 - sinO 


sinhO + sind 


SHELL FORMULAS: The first formula below predicts 
 ІТарве of the shell between frames due to instability. 
The next four predict collapse due to yielding of the 


material according to different yield criteria. 


Formula #92 
von Sanden and Gunther 


Failure due to axial stress exceeding yield stress: 
= 2-06 E а 
Әс D 5 + 1.81К6 


_ Sinhô - sinô 


sinh8 + sin0 


nuoc 
1 + 


To} td 


B and 8 are defined in Formula #88 
Formula 192A 
von Sanden and Gunther 
Failure due to circumferential stress exceeding yield 
stress: 


E 2-8 I 1 
Pc Т E 


Ier HC 
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(1 + M) sinh coss + (1-M) cosh2 sumo 
(ск; E E 





sSinh0 4 sinO 


2 
м= — 3H 


1 - y? 





C is defined in Formula 192 


Wenk's Formula 

This formula was actually developed by Wenk and others 

at DTMB using #92 and #92A to es sas the principle 
stresses. These were applied at mid-bay and mid-thickness 


using the Hencky-von Mises yield criterion. 


h AL 
° a т=т=== 
іс — .5Е + „25 
F = 1 + 20C 


0 0 0 MM, 
sinh; cos; + coshz sins 


sinhO0 + sinô 


Ө апа С are defined in Formutia #92 


Lunchick's Formula 

This formula uses the principle stress relations developed 
by Salerno and Pulos rather than £92 and #92A (9). It 
applies them to the mid-bay location using the Hencky- 

von Mises criterion, however, failure is not predicted 
until the yielding which occurs first in the outer surface 
develops into a fully plastic hinge.  Lunchick's formula 


predicts this yield pressure, but curves (included here 





i pe Е 


as figure A-6) must be used:to determine the actual 


collapse pressure. 


h ЖЕ IA 
Po Бру x: - XZ + 2? 
X = l} + HC 
2 = „5 + 3.62TC 
ge sinh; coss - cosh; sins 


sınh8 + siné 


0, C, and H are defined in Formula %92А 


Entry into the curves requires computation of the 


following ratios: 


Freumterential Bending Stress _ CET 
Circumferential Membrane Stress Im 
Longitudinal Membrane Stress = ЊЕ 
Crrcumferentral Membrane Stress E 


Or = 

4h 
Or = 20g (1 ч 2QC) 
Ong = 20g TC 





SSoJ13g ouegJdquej јетаџодојшподто 
Ssol4g auzaquay [eurpnirsuoT1 | 


SOC UN ттт 
CO г~ xO 


EST 

Sra 
N 

МЖС ЕР. 





EJ 


E de • • ° ° 
==] == = = =] 4 = =i 


Aq одпозодд ртатд 
"за әлпссәлі Ә5481109 





br 
OT 


Stress 


Circumferential Bendin 
Circumferential Membrane Stress 


Curves for Lunchick's Collapse Fressure 
Factor - From DTMB Report 1291 


FIGURE А-6 





gc 


GENERAL INSTABILITY FORMULAS: While Tokugawa's 
formula and graphs of Kendrick's third solution for 
general instability have been given in Appendix A, 
Eryant's simplification of Kendrick's first solution 


has been used in the computer program. 


Bryant's Formula: 





, h а 
Әу 
752 
(n2 — 1 + > ) (п? + А?) 
+ n* - 1 
EI] 
3 
Ra PL 
egies LL > an integer 
by 
n= the integer which gives the lowest per. 


Since out of roundness effects general instability 
collapse, a formula taken from (2) has been used in the 


following form: 





max 








О 


e = eccentricity 


ВА desired collapse pressure 


a pee ond are determined in previous 


formulas. 
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|, EROM PRINCIPLE STRESS CRITERION AT MID-BAY EXTERIOR, 
FORMULA 292A 


1. Formula %92А 


hc 
сү 
Pc ROM e deer 
mm | .85 - B 
1 + @ _ 
2. Y = Pc 
De 
Dm Y 
Y - l Į = 1 
1 + H -85 ~ В 
l + B 
OHNE = É i = 1 
21+ в 


Ha 
Е 
00 
R lu 
T |! 
9 [UJ 
ы 
|| 
о 


H = 0 

a 3 sinh$ coss + cosh% sing 
Sinh + sin 

BEZ au ous d reddas | 


@ = (1.82 L/D_)/\h/Dm 
L/Dg- 2.055 \[h/D,, 


h/D,,= .236(Һ/р,) 7 














3 m i ue | ey 
3 

E 

Ап = 2,055 Dg Sy 

h E 


Equation (A), Section II 





= 75 


Am FROM HENCKY-VON MISES CRITERION AT MID-BAY AND MID- 
THICKNESS, 


VON SANDEN AND GUNTHER FORMULAS 


1. Wenk's Formula 


2ho 1 
БЕ ra 


trj 

| 
= 
+ 
NA 
Ю 

rn 

jo 
Ul 

1 
> 
— 


1 + B 
9 9 9 9 
Q = - sinh 2 cos 2 + cosh 2 sin 2 
sinn O Ft sin O 
2 Y = Pc 
Nn g 
u 
1 
E E 
I  1.1546 Ü. FAL = 
2 4 
2 КЕ ШЕ £ 
F > + 4 = 75 
E Ot ae 
A 1 = 


00 .85 ~ В J 
UC 


ЕЕ = о алоје = а оке а пале 


Ө = (1.82 L/D,) / \\h/D,, 


L/D, = 2.565 |h/Dp 


h/D EZ би) ° 


m 





Equation 





(B), 


Section II 





pe 


SHELL: "IDELY SPACED STIFFENERS 
l. He -von Mises Criterion 
01, , G2 - three principle stresses 
“ 2 2 | 2 2 2 
20.. = (01 05 ) + (05 04 )“ + (0, 01) 
рО Fo: ^ unstiffened tube: -- 
ОТ = 205 05 = 0 
2 E 2 2 2 
20, E 0» +0, E 40, 
B = 305 
G = l = 215076 
2 Т С Í y 
РАЛЕ 6 
y 
m = 2 = 1.1540, 


D o Уе 
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APPENDIX F 


COMPUTER. PROGRAMS 


Two computer programs were used to produce theoretical 
data for this report. Program l utilizes equation (A) 
ОО eeen TEO determine curves of h/D. versus Ау 
for eight values of бу. Since equation (A) is based оп 
the Ө value at which the variable H in formula #92A 
goes to zero, the program first solves for that 0 and 
then gces on to vary h/Dg and Су, solving for Ap. This 
program was used only to define the upper limits of the 
yield failure area and so, is not actually necessary for 
the iterative optimization proceedure of Section IV. A 
flow chart and variable chart are included with the program 


Misting and printout. 


Program 2 follows a scatter-gun philosophy in that 
it determines collapse pressures of given geometries for 
all three modes of failure by various formulas. All of 
the formulas are mentioned in Appendix A and completely 
defined in Appendix E. The input of this program consists 
of outside diameter, desired range of collapse pressures, 
frame proportions, material properties, and range of 
values. Shell thickness is determined by the hoop stress 
relation and frame size by a modified Foppl formula. Program 


2 was not intended to be used to produce a complete 
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optimized design for a given set of requirements. Its 
main use, in solving Eos and printing collapse pressures 
for all modes, is pointing up trends and areas where 
various formulas fail. It has provisions for any 

number of solutions to be printed and any geometry to 

be investigated. Іп short, it is meant to be used as 

a tool by a designer who optimizes, using judgment and 
experience as well as computer time. A flow chart and 
variable chart are included with the program listing 


and a sample of its printout. 





Оооо 


10 


19 


20 


30 
180 


40 


205 
210 


Sea 


FROGRAM 1 


THIS FROGRAM DETERMINES LAMBDA MAX FOR CHANGING H/D 
AND YIELD STRESS USING A FORMULA DERIVED FROM 92A. 


P2 = .09 
NUM = 1 
NEND 
TH = 
DTHP 1 
DTHM ol 
NUM = NUM + 1 

oo. ер) GO TO 30 

Ze 


50 
„7 


H H Kol 


THT = Tå 

DEN = SINH(TH)+SIN(TH) 

EHO = (1.+(3*P2/(1-P2 ))**.5)*S INH(THT)*COS (THT). 
EHI 2 (1.-(3*F2/(1-F2))**.,5)*COSH(THT)*SIN(THT) 


ЕН = -2*(EHO+EH1 )/DEN 
IF(EH) 19,30,9 

NF = NP41 

ТЕ(КР-Е0,3) GO TO 10 

DTHF = DTHP/2 

TH = TH-DTHM 

М = 1 

GO TO 1 

M = Mel 

IF(M.EQ.3) GO TO 20 

DTHM = DTHM/2 

TH = TH + DTHE 

МР = 1 

GO TO ł 

WRITE(6,180) NUM,EH,TH 

FORMAT(' AFTER',13,' ITERATIONS, H IS',F8.6,' AND THETA 
E = 30000000.0 15',26.3) 
DOH = 0,0 

ррон = 40,0 

Крон = 15 

DY = 10000.0 

NY = 8 

DO 50 N = 1.NDOH 

WRITE(6,205) DOH 

Y= 30000.0 

DO 40 M = 1,NY 

ELAM 2 (TH/1.82*DOH*Y/E)**,5 
WRITE(6,210) Y,ELAM 

Y z Y+DY 

DOH = DOH.+ DDOH 

FORMAT(' L4 fBDAS FOR DM/H s ',F7.2,!' ARE ') 
FORVAT(2F12.4) 

END 





EOD 


FRINTOUT OF FROGRAM 1 (Sample) 


AFTER 50 ITERATIONS, H 1S-.000000 AND THETA IS 3.74 


LAMBDAS FOR DM/H - 0.0 АКЕ 
30000.0000 0.0 
40000.0000 0.0 
50000,0000 0.0 
60000.0000 0.0 
70000.0000 0.0 
80000.0000 Qm 
90000.0000 OO 

100000.0000 0.0 

LAMBDAS FOR DM/H = 40.00 ARE 
30000,0000 0.2867 
40000.0000 033311 
50000,0000 0.3702 
60000.0000 0.4055 
70000.0000 0.4380 
80000.0000 0.4683 
90000.0000 0.4967 

100000.0000 0.5235 

LAMBDAS FOR DM/H = 80.00 ARE 
30000.0000 0.4055 
40C00.0000 0.4683 
50000.0000 0.5235 
60000.0000 005735 
70000.0000 0.6194 
80000.0000 0.6622 
90000.0000 0.7024 

100000 .0000 0.7404 

LAMBDAS FOR DM/H = 120.00 ARE 
30000,0000 0.4967 
40000.0000 085735 
50000,0000 0.6412 
60000,0000 0.7024 
70000,0000 0.7587 
80000.0000 0.8110 
90000.0000 0.8602 

100000.0000 0.9068 

LAMBDAS FOR DM/H = 160,00 ARE 
30000.0000 025735 
40000.0000 0.6622 
50000,0000 0.7404 
60000.0000 0.8110 
70000.0000 0.8760 
80000,0000 0,9385 
90000,0000 0,9933 

100000.0000 1.0470 
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A 


EM бӯ 
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( START ) 
| 





pos 0g mEs3. 7 
NUM=1 NEND=50 
DTHP=DTHM=,1 


| NUM=NUM+1 | 
















\ WRITE NUM 


E, DOH, DDOH 
NDOH, DY, NY 





A sr 


EE y = zn 000 | 





= ===>] 


CALCULATE ELAM ELAM 
WRITE Y, ELAM 





L PPD 











TH=TH+DTHP | 


FLOW CHART 


PROGRAM 1 





AN AAA es D TT о панар ч s --- 





moque 


NOMENCLATURE OF PROGRAM 1 








рон - D, /h 

DDOH - Incremental change in Dm/h 

NDOH - Number of changes in Dm/h 

E - Young's Modulus l 

EH - Variable H of yield failure formulas 
ЕҺАМ - x 

P2 - Poisson's Ratio squared 

TH - Variable 9 of yield failure formulas 
Y - Material 6 

py = mu eet change in бу 


МҮ E Number of changes in бу 


EE 





0000000 
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PROGRAM 2 


THIS PROGRAM DETERMINES THEORETICAL COLLAFSE PRESSURES 
DSRESIPIEBFENED CYLINDRICAL PRESSURE VESSELS. THREE MODES 
OF COLLAFSE ARE INVESTIGATED WITH EMFHASIS ON YIELD 
FRECIPITATED BUCKLING BETWEEN THE FRAMES. FRAME SIZE 

1S DETERMINED USING А MODIFIED FOPFL FORMULA. ALL INPUTS 
ARE HANDLED BY MAIN PROGRAM CARDS - NO DATA CARDS ARE 
USED. 


MF = 6 

D = 480.0 

FC 2 300,0 - 
DO 700 NU = 1,14 

Y = 80000,0 

E = 30000000,.0 

Rl = 1,0 

R2 - .11764 

R3 2 1.7 

Р = .3 

F2 2 P**2 

SR = ‚38 

IF(PC.LT.1000.0) SR = ,46 
DSR = .02 

NSR = 20 

T - D*FC/(2*Y) 

DM = D-T 

TOD = T/DM 


WRITE(6,160) T,TOD 

DO 500 IN = 1,NSR 

EL = SR**2*DM*E/Y*(T/DM)**1.5 
FRAMD1 = 12,0 


DFRP = 1.0 
DFRM = 1,0 
МОМ = 1 
NEND = 70 


RATIO = EL*EC/(24*E*MF ) 
KUM = NUM + 1 
IF(NUM.EQ.NEND) GO TO 30 
Е = FRAMDIL*RI] 

B = ЕВАМО] 782 

T2 = B*R3 

W = FRAMDL ~ T2 

YB = YENTRD(B,T,B,W,T2,F) 


YBI = YMOM(B,T,B,W,T2,F, YB) 
DBT = (D-2*(YB+T))%**3 
ВАТ = YBI/DBT 


IF(RAT-RATI)) 19,30,9 


9 NP = NF+l 





10 


1%, 


EO 


30 


ap gc 


ШЕПІР.Е0:3) GO TO 10 
DFRE = DFRF/2 

FRAMD1 = FRAMD1 - DERM 
М = 1 

СО ТО 1 

M = M+1 

IF(M.EQ.3) GO TO 20 
DFRM = DFRM/2 


FRAMDl = FRAMDl + DFRP 

МР = 1 

со TO 1 

YF = YENTRD(0.0 20.0 ,B,4 ‚T2,F) 
Yl = YENTRD(EL,T,B,W,T2,F) 
m s 0, 0, 0,8,4 ТЕСЕ Р" 
Y11 » YMOM(EL,T,B,W,T2,F, Yl 


YBIOYI - YBI/YI 
Y110YI = Y11/Y1 
RONE = (D-2*(T+Y1))/2 
DB = D-2*(T+YB) 
А = F*T2 4 B*"W 
AB = A+B*T 
EB = B*T/AB 
THETA = 10*(12-12*P2)**,25 
TH = THETA*EL/DM/(100*T/DM)**.5 
THT = TH/2 | 
DEN SINH(TH)+SIN(TH) 
= (COSH(TH)-COS(TH))/DEN 
EK = (SINH(TH)-SIN(TH))/DEN 
BETA = 2 ,*EL*EN"*EB/(B*TH) 
EHO = (14(3*F2/(1-F2 ))**.5)*SINH(THT )*COS (THT) 
EH] = (1-(3*F2/(1-EF2 ))**.5)*COSH(THT)*S IN(THT ) 


EH = -2*(EHO+EH1 )/DEN 

EQ 2 -(SINH(THT)*COS(THT )--COSH(THT)*SIN(THT) ) /DEN 
ET - pO Us OHTI-COSR CTHT)SSINCTRT) ) /DEN 
SUBDEN (.85-EB)/(1+BETA) 

EF = l. + 2*EQ*SUBDEN 

EX = 1. +EH*SUBDEN 

EZ = .543.62*ET*SUBDEN 


FTOK = 24*E*Y11/(DM**3*EL*(1-P2)) 

FWIND = 2.6*E*(T/Di4)**2..5/(EL/DM- .45*(T/DM)**.5) 

DENONE = (1+(.85%BETA/EB))/(1+BETA) 

DENUM = 2*T*Y/DM 

FY88 2 26*EX*YBI/(DB**3*BX*DENONE) 

S192 = 1/(.5+1.81*EK*SUBDEN ) 

SI92A - КО вор 

SIWENK = l1./(EF**2-EF/24.25)**,5 

FLUNCH - 2*T*Y/D/(EX**2-EX*EZ4EZ**2)**, 

FBRY s 2*EX*T/DM/850. 4 665E*YlI/(DM^ :J*EL) 

SIGOOR = EC «D*B*DENONE/(AB*2) 

SIGOOR = SIGCOR + E*((D-DB)/2.)*4*(EC/(EBRY-EC) ) /RONE**2 
SIGF = 1000*DM/(4*T) 





500 
700 


902 
903 
904 


922 
923 
924 
225 


1 


E 


SIGT 2 SIGF*2*(1.42*EQ^*SUBDEN ) 

SIGBT = SIGF*2*ET*SUBDEN 

SIGFOT = SIGF/SIGT 

SIGBOT = SIGBT/SIGT 

R= D/2. 

FI = 3.14159 

VT = PI*(EL*(R**2- KBATIALZEBECR-T)552 - (R-T-W)**2 ) 
REX ((R-T -4)**2-(R-T-W-T2)*5*2)) 

WT 2 VT^.284*12/EL 

VOL - EIS (R/12.)**2*]2 

SF92A 2 SI92A*DENUM*VOL/ (WT*100000) 

SFWENK - SIWENK*DENUH*VOL/ (WT*100000) 

Write statements for whatever solutions are required 

SR = SR+DSR - 

РС = РС+100.0 

Format statements 

END 


FUNCTION SUBROUTINES 


FUNCTION YENTRD(EL,T,B,W,T2,F) 

TOF 2 2*F*T2*W + FAXT2X*2' ^ B*W**2 - EL*T?**2 
BOT = 2*(F*T2 4B*W + EL*T) 

YENTRD = ТОР/ВОТ 

RETURN 

END 


FUNCTION YMOM(EL,TLBLWLT2,F,Y) 

FIRST = (EL*T**3 + B%W*%*w3 + FRT2%%*3)/12 
SECOND 2 EL*T*(T/24Y)**2 4 B*W*(Y-W/2)**2 
THIRD 2 F*T2*(T2/24W -Y)**2 

YMOM = FIRST + SECOND + THIRD 

RETURN 

END 





Zol- 
FRINTOUT OF FROGRAM 2 (Sample, continued on next page) 


T IS 3.0INCHES AND T/DM IS 0.00529 
LAMBDA ,EL,FRAME FACTOR 0.380 12.883 4 
SECOND MOMENTS:F,B,1 200.436 439.394 1481.279 
F AREA WINDENBG FORM88 FORM92 FORM92A WENK 
28.2750 -28193.383 590,321 1,857 1.488 1.655 


LAMBDA ,EL,FRAME FACTOR 0.400 14.275 4 
SECOND MOMENTS:F,B,1 225.964 485.456 1680.225 

F AREA WINDENBG FORM88 FORM92 FORM92A WENK 
302022. -42473.102 609.692 1.830 1.469 1.651 


LAMBDA ,EL,FRAME FACTOR 0.420 15.738 4 
SECOND MOMENTS:F,B,1 253.112 533,708  1893.137 

F ARAE WINDENBG FORM88 FORM92 FORM92A WENK 
31.774 90844 ,625 629010425 71,801 1.451 1.637 


LAMBDA ,EL,FRAME FACTOR 0.440 17.273 4 
SECOND MOMENTS:F,B,1 281.881 594,134 2120.124 

F AREA WINDENBG FORM88  FORM92 FORM92A WENK 
33.531 467240.063 647.914 1.770 1.434 1.624 


LAMBDA ,EL,FRAME FACTOR 0.460 18.879 4 
SECOND MOMENTS:F,B,1 312.279 636.729 2351.309 

F AREA WINDENBG FORM88  FORM92 FORM92A WENK 
35.293 62893.883 666.852 1.736 1.420 1.612 


LAMBDA ,EL,FRAME FACTOR 0.480 20.556 4 
SECOND MOMENTS:F,B,1 344.308 691.478 2616.743 

F AREA WINDENBG FORM88 FORM92 FORM92A WENK 
37.058 33034 , 949 685.733 1.700 1.400 1.599 


LAMBDA ,EL,FRAME FACTOR 0.500 22.305 4 
SECOND MOMENTS:F,B,l 377.968 748.369 2886.512 
F AREA WINDENBG FORM88 FORM92 FORM92A WENK 


38.828 22097.688 704.602 1.662 1.325 15867 


LAMBDA ,EL,FRAME FACTOR 0.520 24.125 4 

SECOND MOMENTS:F,B,1 413.264 807.395 3170.714 
FAREA WINDENBG FORM88 FORM92 FORM92A WENK 
40.600 16434.441 723.510 1.622 1.367 1.575 


LAMBDA ,EL,FRAME FACTOR 0.540 26.016 4 
SECOND MOMENTS:F,B,1 450.198 868,543 3468.377 
F AREA WINDENBG FORM88 FORM92 FORM92A WENK 


42.376 pau Tv 782.510 1.581 9950 1563 


LAMBDA ,EL,FRAME FACTOR 0.560 27.979 4 
SECOND MOMENTS,F,B,1 488.769 931.803 3782.55 

F AREA WINDENBG FORM88 FORM92 FORM92A WENK 
44.154 10653.945 761.666 1.538 1.333 1.551 





FLANGE,WEB: LENGTH , THICKNESS 
LUNCHICK RATIOS:F/T,B/T 
LUNCHICK BRYANT OUT ROUND 
1.629 2478.0 54958.6 


FLANGE ‚WEB: LENGTH , THICKNESS 
LUNCHICK RATIOS ,F/T,B/T 

LUNCHICK BRYANT OUT ROUND 
1.610 2526.2 55407.4 


FLANGE ,WEB: LENGTH , THICKNESS 


LUNCHICK RATIOS:F/T,B/T 
LUNCHICK BRYANT OUT ROUND 
1.592 2571.9 55852,6 


FLANGE ,WEB: LENGTH , THICKNESS 


LUNCHICK RATIOS:F/T,B/T 
LUNCHICK BRYANT OUT ROUND 
1.574 2615.4 56292.1 


FLANGE ,WEB: LENGTH , THICKNESS 


LUNCHICK RATIOS: F/T,B/T 
LUNCHICK BRYANT OUT ROUND 
1.556 2656.7 56722 .0 
FLANGE,WEB:LENGTH,THICKNESS 11. 
LUNCHICK RATIOS:F/T,B/T 

LUNCHICK BRYANT OUT ROUND 
1.538 2695.9 57140,1 


FLANGE ,WEB: LENGTH , THICKNESS 


LUNCHICK RATIOS:F/T,B/T 
LUNCHICK BRYANT OUT ROUND 
1.519 2733544 57544 .6 


FLANGE ,WEB: LENGTH , THICKNESS 
LUNCHICK RATIOS:F/T,B/T 


LUNCHICK BRYANT OUT ROUND 
1.501 2769.0 57932,6 
FLANGE,WEB:LENGTH,THICKNESS 12. 
LUNCHICK RATIOS:F/T,B/T 

LUNCHICK BRYANT OUT ROUND 
1.482 2803.1 58301.6 


MANE. WEB Река ТА ,ТАЈСККЕЗ5 


BUSGHNEK RATIOS:F/T,B/T 
LUNCHICK BRYANT OUT ROUND 
1.462 2835.6 58649.1 


9. 


10. 


10. 


10. 


10. 


11. 


PLS 


12. 


80 1.96 7.84 
750 .0080 
WT/FT ЕЕ92 
26118.07 1.246 
10 2.02 8.08 
742 .0095 
WT/FT ЕЕ92 


2565/ ¿38 1.253 


39 2.08 8.31 
7350112 
WT/FT EF92 


792944159 1.257 


68 2.13 8.54 
727 .0131 
WT/FT EF92 


24844 .41 1.265 


95 2.19 8.76 
720 0151 
WT/FT | EF92 
24483.86 1.267 
22 2,25 8,98 
713 01:76: 
WT/FT ЕЕ92 


24149.51 1.268 


49 2.30 9.19 
706 .0198 
WT/FT EF92 


23838.47 1.269 


75 2.35 9.40 
700 .0225 
WT/FT EF92 


23548.48 1.269 


00 2.40 9.60 
604 .0253 
WT/FT EF92 


2900772838. 1.26/ 


25 2.45 9.80 
687 .0284 
WT/FT EF92 


2500235990 1.265 


1T 


EFWENK 
15327 


1.19 


EFWENK 
po 229 


22 


EFWENK 
1.348 


15275 


EFWENK 
1 356 


1:29 


EFWENK 
12203 


1.92 


EFWENK 
1.909 


529 


EFWENK 
1.374 


1.36 


EFWENK 
1.378 


1,81 


EFWENK 
1372 


1.44 


EFWENK 
1.381 





4 


eons 







| MF, D, PC 


a T - 


SR 











| DSR NSR == 
| 


| CALCULATE T, DM, TOD 

\ WRITE Т, TOD 
r= 
| | CALCULATE EL, RATIO | 
| AI 
FRAMEl, DFRP, DFRM, NUM, 
NEND 












! 1A 2A 





YES 


15 
E 3 


| | NO 





| CALCULATE F, B, T2, W, 
YB, YBI, DBT, RAT 


mn YES 
ARA TITO 1 






РЦ 


FLOW CHART 


PROGRAM 2 





zou 


| All calculations used to 
| find solutions to various 
formulas for pe in differ- 


ent modes 









All necessary 


write statements 


| | 


_ END ) 
IS —— 






YES 


= 


IDFRP=DFRP/2 | 


NEU BEN 
Шам к” 
| "=1 | 
| 


lA 


Loop 2 to 2A is the same as 1 to 1A in its method 





— — 


B 

D 

DSR 

E 

EL 

E 
FRAMD1 
MF 

NSR 

P 

PBRY 
ЕС 
РІОМСН 
РТОК 
PWIND 
EXGS 


SF92A 
SFWENK 
SIGFOT 
SICBOT 
2192 
51924 
SIWENK 
SR 

T 

T2 

W 

WT 

VOL 
YENTRD 
YMOM 


e 


NOMENCLATURE OP PROGRAM 2 


Width of frame in contact with shell 
Outside diameter of shell 

Incremental change in slenderness ratio 
Young's Modulus 

Length of shell between frames 

Frame overall width 

Frame overall depth 

ame” Factor 

Number of increments of slenderness ratio 
Poisson's Ratio p 

Collapse pressure due to Bryant's Formula 
Required collapse pressure 

Collapse pressure due to Lunchick 
Collapse pressure due to Tokugawa 5 
Collapse pressure due to Windenburg 
Collapse pressure due to #88 

Frame width/Frame depth 

Web thickness/Frame Depth 

Flange thickness/Web thickness 
Efficiency due to #92A 

Efficiency due to Wenk's Formula 

ср/ст for Lunchick Formula 

opp/op for Lunchick Formula 

u due to #92 

Ú due to #92A 

у due to Wenk 

Slenderness Ratio 

Thickness of shell 

Thickness of flange 

Width of web 

Weight of vessel per foot length 

Volume displaced per foot length 
Subprogram to find neutral axis 
Subprogram to find I 





EO 


APPENDIX G 
COMPUTED DATA 
and 


EXPERIMENTAL DATA FROM THE LITERATURE 





zc 


TABLES OF MAXIMUM VALUES FROM FORMULA #92A AND MAXIMUM 
PRINCIPLE STRESS CRITERION | 


Frame Factor 4 


- 


Frame Factor 5 


Frame Factor 6 


100h/D, 


GJ 
~J 
Oy 


сл 
о 
Uy 


za 


„756 


wol 
2313 
3/6 
‚439 
ОШ 
2:996 
8523 
092 
196 
91.9 
.883 
.946 
1.010 


?) max/1.2 


кҥҥ ҥнҥ ҥнҥ RP Pee 


кҥнҥ = је 


PRR RPP ҥн ыы 


9086 
.9490 
. 26006 
40057 
20253 
.0426 
40570 
«0695 
.0800 
созот 
0987 
.1074 
.1146 


237 


. 9706 
493992 
. 0224 
„0411 
T0970 
-0699 
.0814 
ЕОс 


.9504 
.9864 
¿US 
.0354 
„Оз 
· 05085 
· 0814 
„0316 
. 1016 
„2102 
Pp lis ley a 
.1246 
.1304 


À 


PRR RP RP PRP RP RP PRP Pe 


PRR RPE PRP EPP eee 


PRP RP RPE Pee 


27:16 


52768 


.642 
950 
.485 
.419 


203 
334 


52105 


2276 
ASIS 
.234 
„219 
220 
БОТ 
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TABLES OF MAXIMUM VALUES FROM WENKS FORMULA 


ee ж-- -æ 


HENCKY-VON MISES YIELD CRITERINN-NIDPLANE, MIDBAY 





100h/D,, 7 пах/1.2 A ү, 
Frame Factor 4 4 2925 251333 222 16426 
„439 1.164 · 76 Iso 

503 1.198 272 1562 

‚566 TRIP КӨТ 128523 

‚629 2 7 63 1.506 

„692 1.248 ‚60 1.481 

. 756 T. 265 753 1.454 

. 819 47270 455 1,438 

ES 1.280 52 12425 

.946 22292 50 1.409 

1.010 192305 ‚49 J. 391 

Frame Factor 5 3716 1.1592 84 196620 
.439 1212390 76 L 2571 

- 503 1.2139 270 12,538 

2566 10523247 .67 1.497 

2629 1.2500 263 1.472 

‚692 1.2644 ‚60 1.448 

1756 1.2759 2057 19428 

Frame Factor 6 23.26 1.1794 ‚84 573 
‚439 1.2068 276 1.540 

.503 1.2298 . ZL 15502 

. 566 ІІ 247 .66 1.476 

.629 1.2609 .62 [m 

.692 1.2744 ‚60 ПИА АЗ 

Эб 1.2859 51 1405 

.819 1.2960 254 1.390 

.883 1,3032 15 1553774 

.946 15311983 DEL 170350 

1.010 1, 34076 ‚49 1,345 





FRAME PARAMETERS AT 


-09- 


elu EFFICIENCY 


%92А 


EODD. 1]/Ip 100b/D 
Frame Factor 4 „251 80 4.758 2105 
513 272 52555 2350 
316 ‚ 66 5.860 ‚2485 
· 439 .61 6.340 ‚2640 
„503 Eg 6.817 27120 
966 24 7 3238 .2780 
.629 Sil: 1.6055 22050 
.692 49 82072 .2930 
эб 46 8.464 - 2960 
. 819 ‚45 BB ‚3045 
. 883 .44 9.267 . 3120 
. 946 .42 9,662 . 31.45 
1010 40 1702065 ‚3150 
Frame Factor 5 „251 ‚80 5.064 ‚2065 
9313 “72 52672 „2215 
‚376 . 66 6.224 „2340 
· 439 .61 6.735 .2420 
‚503 . 57 7 2212 . 2525 
‚566 uo 7.679 ‚2610 
2029 51 8.124 20675 
‚692 ‚49 So. 567 22005 
1756 КО] 9.001 ‚2810 
Frame Factor 6 5 Т 82 523564 . 1990 
313 E72 5.958 ЛІМІ 
72715 .66 6:539 ‚2220 
439 .61 7.062 2315 
2503 234 7.566 .2400 
‚566 254 8.056 22475 
. 629 51 825287 2535 
‚692 ‚49 8.995 22:505 
. 756 ‚47 9,451 ‚2660 
‚819 Nl 9.919 52700 
.883 243 10.391 2207125 
.946 222 E0853 „2785 
1.010 ‚40 11,354 ‚2790 
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FRAME PARAMETERS AT MAXIMUM EFFICIENCY #WENK 
100h/D, Iy/Ir 100b/D 
Frame Factor 4 276 292. 6,292 ‚2800 
439 27b 6.647 ¿2925 
503 272 72085 3055 
566 67 7.468 > 140 
629 63 7.8413 . 32106 
692 ‚60 B2] ‚3290 
27:36 238 5.580 22915 
819 55 8.930 . 3420 
.883 52 9.280 . 3440 
. 946 оО 9.634 „3495 
IE .49 9.986 «9 95/5 
Бете Factor 5 2570 · 84 65068 52670 
. 439 . 76 7.008 (2755 
503 . 70 7.2453 ‚2830 
„566 67 715678 ‚2950 
‚629 263 8.264 ‚3010 
‚692 ‚60 85657 .3090 
756 а 9,041 SOBALD 
Frame Factor 6 a36 .84 67555 -2540 
.439 76 TRS E „2620 
E503 ТШ TEES 2715 
. 566 .66 8.195 ‚2780 
‚629 HO» 8.616 22855 
.692 ‚ 60 9,087 „2930 
2471512 КЭ 9.446 .2980 
.819 ИЕА 9.859 . 3010 
. 883 452 T2 9 .3060 
.946 Кој 70.678 +3130 
OTO . 49 12.21.06 - 3160 
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EXPERIMENTAL DATA FROM REFERENCE 9 


Model Number  100h/D Jy/1000 Х Woon Wexp 
1 .88 46.8 .41 1.43 1.68 

2 „784 47.5 „42 1.42 je 

BR7m . 784 59.2 ‚46 1.44 1.62 
BR7 „599 44.8 „59 ШОТТ dS 

„814 47.1  - .61 ol 1.24 

„474 49.6 67 1.09 1.24 

0 70 .997 ШТІ! 


7 5291 47. 


W exp represents Y determined using actual collapse 
pressure and Woon was determined using collapse pressure 


calculated by formula 792A. 
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EXPERIMENTAL DATA FROM REFERENCE 19 


General Model Parameters 


.1405 


I 


Left/D, 


E 30,000,000 ib/in? 


| 


Light Frames 


MODEL NUMBER 100h/D  Óy/1000 ^ V exp 
79 1301 27 .876 „92 
82 „303 28 „855 „90 
76 „307 42 1.06 ‚ 66 
88 „293 27 „890 „85 
75 „294 39 1.06 „91 
81 „295 28 „900 „93 
84 „300 29 „909 .81 
87 .300 29 .905 УБ 
80 0 28 .840 .87 
85 2:31] 29 .884 .83 
92 .288 32 1.00 .85 
91 .290 32 .975 .84 
93 .308 29 .886 .88 
96 .298 28 .896 .85 
94 2298 28 . 896 . 84 
90 .300 27 .88 „82 


In the computation of ^ values, h/D was used instead 
p D resulting in a small error which was not critical 


in the graphical methods used. 


ez represents the ratio of actual collapse pressure to 


the hoop stress collapse pressure. 
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ALUMINAUT 


in the following section, A, y, and n have 
been found for the pressure hull in Aluminaut. The 
dimensions of that craft were taken from References 


(5) and (15) and are for the cylindrical hull alone. 


Operating Depth = 15,000 ft. | LOA = 33.3 ft. 
Safety Factor = 1,4 nc - OLD EET 
oy = 60,000 psi D = 8.0 ft. 
E = 29,000,000 psi ры = Sort, 
p Eom h = Str 
NT _ | ЕКЕ ЖТ 





(h/D,,) * E 


А = .324 





- 1.4 Х 15,000 Х .446 
9400 psi 


~ 
= 
|| 
rd 
Q 
| 


^59y | Pc 





Ww 10' 


ү = 2.65 1n. 


From the above computations, it is apparent that the 
designers of Aluminaut produced a vessel whose parameters 
put it in the low lambda design region, However, it has 


been noted (15) that a more efficient structure is possible. 
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